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Title : Method for rdentification of Algorithm for Analyle Delerminaiio^ 

FTELD OF THE INVENTION 

The present invention relates to the fields of spectroscopy^ 
spectrophotometry^ and chcmometrics. In particular, the present 
5 invention relates to a method of calibrating a spectroscopic device for 
predicting analytc levels in a non-invasive manner, llie metliod of the 
present invention is particularly suitable for blood glucose prediction 
based on near^infrared spectrophotometry measurements. 

BACKGROUND OP THE INVENTION 
10 Biotechnological analysis and examination are oftcm based on 

the measurement of various chemical analytcs in the composition of a 



{3i biological matrix such as blood, interstitial fluid, or living tissue* Such 

measurements may be used to evaluate a patient's state of health and to 
f J determine what, if any, treatment is necessary. For example, the frequent 

iiJ? 15 monitoring of blood glucose levels in diabetic persons with glucometers is 

often necessary to allov^ such persons to manage the diabetes mellitus 



disease, by taking insulin injections or oral drugs to lower blood glucose 



LJf when required. Intensive treatment based on frequent blood glucose 

il. 

^j.; measurements can significantly reduce the incidence of blindness^ kidney 

h3c 20 loss/ and other complications associated with diabetes. 

Most home-based glucose measurement systems require the 
patient to invasively collect a blood sample, by pricking his or her finger, 
placing the sample on an appropriate test strip, and then testing the 
sample in an optical glucumeter. For millions of diabetics around the 
25 world/ the use of lancets or other sharp instruments to draw blood for 
monitonng their insulin levels is a painful process, and one that often 
builds up calluses on fingers, making the collection of blood even more 
difficult. This Invasive procedure may be especially difficult to perform on 
children and therefore particularly trying on parents. Furlhermore, the test 
30 strips required for each blood sample are generally not reusable, and when 
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multiple measurements are taken each day, amount to significant costs 
from the patient's point of view. Thus, despite the fact that a large number 
of diabetics should take several measurements Ihroughoul each day (for 
some individuals, physicians reaimmend testing glucose levels from 4 to 

5 7 limes daily), due to the pain, cost, and inconvenience involved, many 

dinbcticp do not monitor their glucose levels freq uently enough. 

Spectroscopy is based on the analysis of how incident 
radiation interacts with the vibrational and rotational states of molecules 
which arc of analytical interest. Spectroscopic measurement techniques 

10 have gained increased fK>pularity because of the ability to provide fast and 
non-invasive measurements of concentratioiw of different chemicals or 
analytes. For the reasons indicated above, this is particularly desirable for 
home based glucometers. Speclropholomelry is a type of spectroscopy 
commonly used to quantitatively measure concentrations of analytes 

15 based on spectral energy distribution in the absorption spectrum of a 
sample solution or medium. In spectrophotometry, the energy 
distribution is typically analyzed within a range of the visible, ultraviolet, 
infrared, or near-infrared spectra. For example, near-infrarcd radiation 
(NIR) is electromagnetic radiation having a wavelength of between about 

20 0.75 and 2.5 micrometers (i.e. from 150 to 400 THz), Near-infrared 
Spectrophotometry generally uses instruments with quartz prisms in 
monochromators and with lead sulfide photoconductor cells as detectors 
to observe absorption bands, and NIR spectrophotometry is increasingly 
being used to measure in vivo analytes such as glucose, fructose, glycerol, 

25 and etfumol. 

Non-invasive, spectrophotomctric measurement of glucose 
in human beings is performed by focusing an incident radiation source (or 
sources) on a specific part of the body and detecting the spectral 
distribution of the radiation transmitted therethrough. The absorbancc of 
30 light from the incident radiation is due to the chemical components 
within that body pari including water, fat, protein, hemoglobin, melanin, 
and glucose. One difficulty with gluaise measurement spectral analysis, is 
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the spectral overlap between glucose and other chemicals foimd in blood, 
often in much greater quantities than glucose. In addition, the thickness, 
color, and structure of the fikin, bonds^ and blood through which the 
incident radiation passes will affect the transmission (or 
5 reflection/absorplion) thereof. Furthermore, the concentration of analytcs 
can vary with changes in activity level, diet, hormone fluctuations, and 
other factors. Glucose concentration measurements arc also particularly 
susceptible to variations in physical and chemical conditions including 
temperature, pressure, humidity, and skin hydration. As a result, to 

10 perform a reliable non-invasive glucose prediction, NIR spectral 
measurements should be performed through a vascular equilibrated 
rcgion of the body, and a NIR glucose spectrophotometer must be carefully 
designed so that the quality of raw spectral information from an NIR 
glucomctcr is high. See generally Waynant and Chenault, "Overview of 

15 Non-Invasive Optical Glucose Monitoring Techniques", IEEE LEOS 
Newsletter, voL 12, no. 2 (April 1998); and Burmcistcr and Arnold, 
"Spectroscopic Considerations for Noninvasive Blood Glucose 
Mcasutxjmcnts with Near Infrared Spectroscopy", IEEE LEOS Newsletter, 
vol. 12, no. 2 (April 1998). 

20 Near-infrared glucose measurements arc generally suitable 

for tissue depths ranging from about 1 mm to 10 cm, and arc often 
performed through a patient's finger tip, although other areas of the body 
(for example the web between two fingers, an car Jobc, or the upper lip) can 
also be used. The sample thickness is an important experimental 

25 parameter because a greater thickness increases the amount of absorption 
and thereby lowers the minimum limits for detection, whereas because 
less incident radiation successfully traverses through a thicker sample (i.e. 
without being absorbed) effectively increasing the spectral noise: see 
Burmcistcr and Arnold, supra. 

30 In addition of the difficulties in obtaining accurate spectral 

measurements with NIR spectrophotometry, a more significant difficulty 
associated with such measurements has been the need to calibrate such an 



instioinicnt so that it may be used by various different individuals, whose 
analytc concentrations and variability^ for example, may differ 
substantially. Calibration of spectrophotometers (and of analytical 
instruments in general) is necessary to ensure the accuracy of 

5 measurements performed by such devices. 

Two approaches to calibrate and predict anal yte (e.g> glu cose) 
levels have been used in the prior art. In the first approach, a thorough 
calibration process is performed separately on each individual with whom 
the instrument is to be used. This individual calibration process requires 

10 taking a number of blood samples over a period of time from the 
individual^ and obtaining reference glucose concentration measurements 
from these samples, A number of corresponding non-invasive 
spectroscopic measurements are taken concurrently, and calibration 
regression analysis is then performed to correlate, e.g. through linear 

15 regression analysis, the non-invasive measurements with the reference 
measurements. In this manner, the spectroscopic device or 
spectrophotometer is custom calibrated to the specific individual. A major 
disadvantage of this approach, however, is that the calibration model is 
limited to data from that particular individual which may be an 

20 insufficient amount of data for highly accurate analj^e prcdictiom. To 
alleviate this problem, the custom calibration process may have to be 
carried out frequently, even on a daily basis. 

An alternative calibration approach requires that a universal 
or general calibration algorithm be developed. This approach is based on 

25 the notion that by making a large number of calibration measurements (as 
in the first approach) a single calibration model that sufficiently accounts 
for all individual variability can be calculated. Chemometrics, the 
application of mathematical, statistical and formal logic methods to 
chemistry^ is generally used to process and compute the spectral intensity 

30 data and to produce a calibration model: see Small and Arnold, ''Data 
Handling Issues for Near-Infrared Glucose Measurements'', ILEE LEOS 
Newsletter, vol. 12, no. 2 (April 1998); and Shaffer, Small, and Arnold, 



"Genetic Algorithm -Based Protocol for Coupling Digilal Filtering and 
Partial Uast-Squarcs Regression: Application to the Near-Infrared 
AnalysiB of Glucose in Biological Matrices," Anal Chem. 68, 2663-2675 
(1996). However, for measurements of glucose, and other analyies which 
vary distinctly and differently from one individual sample to another, a 
single universal calibration algorithm is often ineffective and can result in 



significant, and sometimes dangerous, erroneous predictions of analytc 

concentration levels. 

Thus, while NIR measurements provide a non-invasive, fast, 

10 painless, and convenient technique to monitor glucose levels, correlation 
and clinical interpretation of spectral measurwnents to obtain the true 
glucose levels is crucial for proper therapy and disease management. 
Proper calibration of instruments for different patient populations (which 
wUl vary in ethnicity, age, weight, and so on) is crucial in obtaining 

15 accurate glucose prediction models. Furthermore, careful validaHon and 
testing of the non-invaslvc results and the glucose prediction equation is 
needed to determine if the glucose correlation is consistent in all clinically 
important conditions and for aU, or at least most, types of patients. 

fiyi yiMARY OF TH F IMVFNTION 

20 In a first aspect the present invention provides a method of 

calibrating a spectroscopic device for providing a non-invasive 
measurement of an analytc level in a sample, comprising: (a) generating a 
plurality of calibration algorithms; (b) talcing a set of non-invasive 
measurements on said sample with said spectroscopic device; (c) 

25 calculating a predicted set of analyte levels for each of the caUbration 
algorithms in i«sponse to the set of non-invasive measurements, each of 
tfic predicted sets of analyte levels being characterized by a slope, a 
varlabUity range, and a standard error of prediction; and (d) selecting an 
appropriate calibration algorithm based on the slope, the variability range, 
30 and the standard error of prediction for each of the predicted scte of analyte 
levels. 



Preferably^ step (d) comprises (i) selecting the predicted sets of 
analyte levels in which the standard error of prediction is less than an 
upper error limit, the variability range is greater than a lower range limil/ 
and the slope is between a first lower slope limit which is less than 1 and a 
first upper slope limit which is greater lhan 1; (ii) for each of the predicted 
sets selected in step (i), calculating a suitability score in response to th e 
slope and the istandard error of prediction for that predicted set, and 
selecting the calibration algorithm corresponding to the predicted set 
having the optimal suitability score as the appropriate calibration 
algorithm; (iii) if no predicted sets arc selected in step (i), selecting the 
prcdictoj sets of analyte levels in which the standard error of prediction is 
less than the upper error limit, the variability range is greater than the 
lower range limit, and the slope is between a second lower slope limit and 
a second upper slope limit which is less than or equal to said first lower 
slope limit; (iv) for each of the predicted sets selected in step (iii), 
calculating a suitability score in response to the slope and the standard 
error of prediction for that predicted set, and selecting the calibration 
algorithm corresponding to the predicted set having the optimal suitability 
score as the appropriate calibration algorithm; (v) if no predicted sets are 
selected in step (i) or step (iii), selecting the predicted sets of analyte levels 
in which the variability range is lower than the lower range limit and in 
which the standard error of prediction is less than the upper error limit; 

(vi) from each of the predicted sets selected in step (v), selecting the 
calibration algoritfim conesponding to the predicted set having the lowest 
standard error of prediction as the appropriate calibration algorithm; and 

(vii) if no predicted sets are selected in step (i), step (Ui), or step (v), 
determining thai no calibration algorithm is appropriate. 

Also preferably, step (a) comprises (i) compiling non-invasive 
and corresponding reference data sets of measurements of analyte levels 
for each of a number of samples; (ii) disregarding data sets that are not 
suitable for calibration; (iii) combining data sets that are suitable for 
calibration into a plurality of groups depending on whether correlations of 
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ihe combined data sets meet predetermined criteria; and (iv) generating a 

calibration algorithm for each of the groups of data sets. The criteria in step 

(iii) may inchide minimizing or lowering correlations of the combined 

data sels in a parlicular group with parameters other than said analytc and 
5 maximizing or increasing the correlation between data sets in a partiadar 
group. Convenienily, steps (iii ) and (iv) are performed using partial least- 

squares regression analysis. 

Sleps (a), (c)^ and (d) of the method of the present invoition 

may be performed on a computer associated with the spectroscopic device. 
10 FurlliermoK?, the analytc may be glucose, the sample may be an individual 

patient^ and the spectroscopic device may be a near- infrared 

spectrophotometer. 

The objects and advantages of the present invention will be 

better understood and more readily apparent with reference to the 
15 remainder of the description in conjunction with the accompanying 

drawings. 

ff RTEF PESCRlPnO N CT? THE DRAWINGS 

In the drawings which illustrate, by way of example, 
preferred embodiment of the invention: 
20 Figure 1 is a correlation scatterplot of glucose predictions for a 

number of different patients based on a universal, single calibration 
model; 

Figure 2 is a correlation scatterplot of glucose predictions for a 
number of patients based on the multiple algorithm calibration model 
25 according lo the present invention; 

Figure 3 is another correlation scatterplot of glucose 
predictions for a number of patients based on the multiple algorithm 
calibration model according to the present invention; 

Figure 4 is a correlation scatterplot of glucose predictions for a 
30 typical patient based on a tmiversal, single calibration model; and 

Figure 5 is a correlation scatterplot of glucose predictions for a 



typical patient based on the multiple algorithm calibration model 
according to the present invention. 



The present invention can be used with a typical NIK 
_5.. spcvctto photomctcr system having a light source which is projecte d 
through the item to be examined, a sample interface mechanism, a 
spectrometer to separate the light into its component wavelcngdis, a 
detector, amplification electronics and a computer. By measuring the loes 
(absorption), between the source and the detector and applying appropriate 
10 chemometric (mathematical) techniques, it is possible to non-invasively 
determine the chemicals being examined since different chemicals absorb 
different amounts of light. Such a spectrophotometric device and method 
are described in detail in United States Patent No. 5,361,758. the contents of 
which are incorporated herein by virtue of this reference. While the 
15 present description relates primarily to glucose measurement, one of the 
major fields of application for NIR measurement at present, it will be 
understood that the principles of the present invention cquivalently apply 
to other analytes non-invasively measured using various spectroscopic 
techniques. 

20 To utilize the NIR spectrum for glucose measurement, it is 

necessary to use a spectrometer which has wide dynamic range, a high 
signal to noise ratio, and exhibits low scattering losses. The output from 
the spectrometer is used to generate spectra with high precision both in 
absorbancc and wavelength. Significant glucose absorption bands are 

25 centered about wavelengths of 1-67, 2.13, , 2.27, and 2.33 micrometers (as 
discussed in SmaU and Arnold, "Data handling Issues for Near-Infrared 
Glucose Measurements", supra). To be able to use NIR to measure a 
particular compound/analyte, chemometric mathematical analysis is 
applied to the measured spectrum. The mathematical analysis techniques 

30 are carried out by a ajmputer equipped with advanced software capable of 
interpreting the resulting complex spectra. 
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To uiiivcrsally calibrate an NIR spectrophotometer^ Ihe 
absorbance associated with the compound or analyte of interest mwt first 
be measured on a relatively large number of samples. These NIR 
measurements arc then compared to measurements made in a more 
5 traditional and more accurate invasive manner. From lhe«Je comparisons 

an algorithm is develo ped that characterizes the analyte to be measure d. 

The methods used to generate calibration models for NIR 
spectrophotometers arc sophisticated. An important criteria is that 
calibration samples be uncorrelated with other chemical species 
10 concentrations- Satisfaction of this criteria is burdensome because of the 
large number of chemical species typically present in a biological matrix or 
growth medium. 

m The present invention provides a multi-algorithm 

l^l calibralion method for use with a spectrosci^ic device capable of predicting 

J J 15 analjrtc levels in a biological matrix, such as blood glucose levels, for a 

l5t variety of different samples with significantly greater accuracy tf\an prior 

art universal single calibration models. The calibration method of the 
invention has two key parts: the generation of a plurality of different 
l^i algorithms and the selection of the algorithm which is most appropriate 

1^' 20 for gmerating an accurate prediction of analyte concentration. 

To generate suitable calibration algorithms, the non-invasive 
and reference measurement data sets for a number of individuals are first 
compiled. Once this is achieved, correlation techniques are carried out to 
relate non-invasive spectral measurements to reference glucose levels. At 
25 this stage the suitability of individual data samples for calibration purposes 
must be considered. If some of the non-invasive spectral measurements in 
a data set result in a saturation of the measurement equipment, e.g. at the 
output of an analog-to-digital converter, then the data sample should be 
disqualified for use in a calibration algorithm. Similarly, if the range of 
30 variability of the non*invasive spectral measurements in a data set is too 
small, the data sample should also be disqualifxed from calibration. 
Furthermore, if the correlation of glucose in a particular data set is 
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undesirably correlaled to another variable, such as the concentration of 
another chemical, the sample should again be disqualified. As a result, a 
large numbers of sample sets will not be suitable for calibration algorithm 
development. 

5 The grouping or combining of calibration-suitable data sets 

{each gpocific to a particular in dividual) into different al gor ithms oc curs 
on a trial and error basis, by repeatedly combining these data sets and 
performing a new correlation on the combined sets of data. If the 
correlations of the combined sets meet certain criteria, then the 

10 combination is acceptable, if not a different combination is attempted. 
Once a combination is deemed acceptable, the addition of a further data set 
is attempted and its acceptability is determined in turn. In this marmer, a 
combination of data sets for an algorithm grows from 2 data sets, to 3 data 
sets, to 4 data sets, and so on until a sufficiently large number of data sets is 

15 obtained from which a reliable calibration algorithm can be developed. 

In combining data sobs for a single algorithm, two criteria arc 
key. First, any correlation of glucose in the combined data set group with 
other parameters must be minimized. Such non-desirable correlations 
may be exhibited with different chemical species concentrations or other 

20 time varying quantities. A suitable quantitative test for this criteria is to 
ensure that all undesirable correlations of the combined data sets, as 
indicated by the square of the correction coefficient are less than a 
certain limit, such as R2 < 0.25. In known manner, the correlation 
coefficient, R, is a measure of the tendency of two variables to vary 

25 together. If the addition of a data set to another data sot or to a suitable sub- 
combination of data sets docs not satisfy this criteria, the new combination 
will not be acceptable- 

In addition, the combination of data sets can be optimized by 
ensuring that the data sets exhibit a high degree of correlation with one 

30 another, so a second criteria in combining data sets should require that the 
data sets in a particular group be sufficiently correlated with one another. 
This can be achieved, for example, by ensuring that the square of the 
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cori^elatjon coefficient is greater than about 0.9 for an existing sub- 
combination of data sets and another data set being potentially combined 
with that sub-combination. 

In an alternate embodiment of the present invention, in:>tead 
of succcsBivcly adding on to, or building up, groups of data sets for a 
ca libra tion algorithm, a large group of calibration suitable da ta sets can be 
initially taken up at random. In this embodiment, data sets are removed 
from the group, and the correlation criteria above arc reevaluated. If the 
correlation criteria improve after the removal of the data set the data set is 
permanently removed from the groups whereas if the criteria do not 
improve, the data set is reinstated back into the group. Using this top- 
down approach, the removal of data sets can cease when the critma meet 
acceptable limits. 

The correlation calculations for determination of the 
calibration algorithm groupings (and ttie associated calibration algorithms) 
become increasingly complex with increasing amounts of data and arc 
preferably carried out using partial lcast«squares (PIS) regression analysis. 
The PLS analysis techniques arc well known to those skilled in the art and 
provide good linear approximation by removing considerable redimdant 
information: see generally Small and Arnold, "Data handling Issues for 
Near-Infrared Glucose Measurements", supra. Other types of analysis, such 
as partial component analysis (PCA) or artificial neural networks (ANN), 
can also be used. 

This calibration algorithm generation process is used to 
generate at least two groups of data sets (corresponding to groups of 
individuals), each providing a calibration model or algorithm for 
predicting the level of analyte present in a biological matrix based on a 
non-invasive spectroscopic measureimsnt. Preferably, two initial groups of 
data set arc developed^ and, from these two groups, additional groups can 
be generated. These additional groups are generally smaller in size then 
the two initial groups and may meet more stringent cocrclaticm criteria. 
These additional groups may be constructed with the top-down approach 
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described above, wherein data sets arc sclccUvcly reinoved from an initial 
group to determine if the correlation crileria improve or deteriorate after 
the removal. 

Given the availability of multiple calibration algorithms, it is 
5 necessary to select the most suitable of these algorithms for predicting the 
levels of analyte in an y g iven patient or sample. O n ce this is done t he 
spectroscopic device is calibrated with the selected algorithm. The 
algorithm selection process is generally independent of the process used to 
generate the multiple algorithms, although it will clearly depend on the 
10 number of different algoritiims that are available. 

l»or any given patient or sample on which the non-invasive 
spectroscopic measurement is to be made, it is necessary to firsl compile a 
set of highly accurate reference measurements and corresponding non- 
invasive spectral measurements at selected intervals over a period of time. 
15 In a preferred embodiment, two successive non-invasive spectral 
measmements arc successively taken each measurement interval, and 
their mean is then taken as a single non-invasive measurement for that 
interval. 

Before algorithm selection begins, two variables may be set to 
20 determine the stringency of acceptance requirements for patients (i.e. 
whe^er at least one of the avaUable calibration algoritiims will be found 
suitable for a patient). These are (i) a lower range limit for the range of 
glucose or analyte variability in the data set (if the data does not vary 
signiHcantly it may not provide sufficient information to be useful), and 
25 (ii) an upper error limit (SEP„„) for the standard error of prediction (SEP), 
i.e, the square root of the mean of the squared deviations of the reference 
analyte concenh^tions from the analyte concentration values predicted by 
a calibration algorithm. 

From the non-invasive spectral measurements, a predicted 
30 set of analyte concentration levels are calculated according to each 
calibration algorithm. When compared to the reference set of 
concentraHon levels, each predicted set will be calculated from a linear 



-13' 



regression or "besl fit" line characterized by: a range of glucose level 
variability; a slope (which ideally equals 1, i.e. the predicted values 
correspond directly to the reference values); a correlation coefficient 
defining the correlation of the predicted value set with the reference value 
5 set; and a standard error of prediction (SEP), as defined above. Note that as 
the cor relation coefficient R increases the SEP will decrease, and vice ver sa. 

The algorithm selection process can include 4 levels or steps. 
In a first level, if the predicted data set for an algorithm meets the 
following criteria 

i) SEP < upper error limit (SEPmax) 

ii) Analytc range > lower range limit 
and iii) 03 < slope (of regression line) < 1.05, 

a suitability score is calculated according to tfie foUovdng equation: 

score = (slope) (R2) / (SEP) 

If more than one algorithm meets the first level criteria, the algorithm 
with the highest suitability score is selected and the selection process is 
complete. If no algorithm qualifies in terms of the first level criteria, the 
selection process advances to a second level. 

The second level is identical to the first level, except that the 

20 slope criteria becomes 

033 < slope (of regression lirc) < 05 

For algorithms that meet this revised slope criterion (and the SEP and 
range criteria), a suilabiUty score is calculated as in the first level. Again, If 
more than one algorithm quaUfies at the second level, the algorithm with 
25 the highest score is selected, and the selection process is complete. 

At a third level, if none of the algorithms qualify at either the 
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first level or the second level and the range of analyte measurements is 
greater than the lower range limit, the sample or patient is excluded from 
mcajsurcmtint^ :>ince under these circumstances no calibration algorithm 
can be relied on to perform accurately. Lastly, at a fourth level, if none of 
5 the algorithms qualify at either the first level or the second level and the 
range of analyte measu re ments is lower than the lower range limit, the 
algorithm with the lowest SEP is selected if that SEP is less than the upper 
error limit. If the SHP is greater than the upper error limit, the patient is 
again excluded from measurement. 
10 Note that the slope criteria limits in the first and second 

levels described above are preferred limits, and these levels can also be 
varied somewhat (similar to the lower range limit and upper error limit) 
^ff, depending on the application and acceptance requirements. 

To illustrate the principals of the present invention, the 
tH: 15 multiple algorithm calibration method was used to calibrate a NIR 

f !^ spectrophotometer for measuring glucose levels in a number of different 

^ palienls. All reference data measurements were taken with the industry 

standard YSI (Yellow Springs Instruments) Glucose Analyzer which 
provides precise glucose measurements on (invasive) blood samples. The 
idl 20 inventor found that only approximately 6% of NIR patient data sets 

qualified for calibration purposes, given the criteria described above. These 
calibration-<iualified sets of data were grouped into two initial calibration 
algorithm sets, each initial group was built up to a group with 33 sets of 
data within it. Eleven additional groups were sptm off from these two 
25 initial groups and each of rticsc groups had at least 8 sets of data within it 
Thus, thirteen different calibration algorithms were available for 
pivdicting glucose levels. 

The algorithm selection was based on a total of 24 reference 
(i.e. YSI measured) and 24 means of two NIR finger measiurements, taken 
30 from a palient over a three day period. The above described algorithm 
selection process was run for a number of different patients (not involved 
in the calibration process). 



-15- 



Figure 1 shows a corrclatjon scatterplot of glucose predictions 
for a number of different patients based on o universal calibration model 
(QUAD 37). The scatterplot in Figure 1 is superimposed on a Clarke error 
grid, as described in Clarke and Cox, "Error Grid Analysis'^ Diabetes Care, 
5 10:622-628 (1987). The Clarke error grid breaks the correlation space into 
five regions (A-E) that assess measurement accuracy on the basis of 
validity of the corre^sponding clinical decision. Correlation points falling 
within the "A" region correspond lo the correct clinical decision being 
made based on the similarity between the actual and predicted glucose 
10 levels, (i.e. tiie predicted values deviate by no more lhan 20% from the 
reference values). In the ^B"" region, the predicted values deviate by more 
than 20% from the reference values^ but treatment decisions made based 
on the predicted levels of glucose would not jeopardi/x! or adversely affect 
the patient- In regions "C", "D", and "H", the predictions significantly 
15 deviate from the reference values, and treatment decisions based on these 
predictions may well be harmful to a patient The Clarke error grid is often 
used to evaluate the clinical consequences of home-use blood glucose 
monitor errors, in contrast to more conventional analytical methods 
which may be more suitable for laboratory reference devices. 
20 Referring to Figure 1, it can be seen that the tipper error limit 

SRPmax universal calibration algorithm is 3.07 mmol/L and ttvat 

significant amounts of the data are outside acceptable limits or regions. In 
contrast/ Figure 2 shows the correlation scatterplot of glucose predictions 
for 25 patients* based on the multiple algorithm calibration model 
25 according to the present invention, with SHP^a^ equal to 2.76 mmol/L. 
Similarly, Hgurc 3 shows a correlation scatterplot of glucose predictions 
for 21 patients based on the multiple algorithm calibration model 
according to the present invention, with SEP^ax equal to 2.28 mmol/L. 
Both the predictions in Figure 2 and particularly Figure 3 exhibit much 
30 better accuracy than those in Figure 1 (calibrated according to a universal 
algorithm). 

Similarly, Figure 4 shows a set of correlation data for a typical 
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patient, calibrated according to the same universal algorithm, whereas 
Figure 5 shows the correlation data for the identical patient calibrated 
according to the multiple algorithm method of the present invention. 
Once more, the impmvcd accuracy of the glucose prcdictioixs based on the 

5 calibration method of tlie present invention is evident. 

It will be clear to those skilled m the art that the principles of 
the present invention arc applicable to a broad range of spectroscopic 
applicatioiu;. For example, NIR radiation can be used to penetrate such 
items as human tissue, vials of blood, or containers of milk - all normally 

10 considered opaque. Mathematical analysis of the resulting absorption 
spectrum determines the composition of the substance penetrated by the 
light. Thus, in a blood testing lab, samples can be examined without 
contact by passing the light through a plastic vial containing the liquid. 
Similarly, in a dairy application, the NIR Ught is passed through a bag of 

15 milk and such parameters as butlerfat, solid particles and lactate arc 
measured. The multiple algorithm method of the present invention is 
easily extendible to these types of non-invasive measurements as well. 
Furthermore, the method of the present invention can also be applied to 
the non invasive measurement of a number of other blood analyles such 

20 as cholesterol, hemoglobin, HblAc, Frudosaminc, and 1.5 AG. 

While preferred embodiments of the present invention have 
been described, the embodiments disclosed are illustrative and not 
restrictive, and the invention is intended to be defined by the appended 
claims. 
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I CLAIM: 



1. A method of calibrating a spectroscopic device for providing a 

non-invasive measurement of an analyte level in a {sample, comprising: 
(a) generating a plurality of calibration algorithms; 
5 (b) taking a set of non-invasive measurements on said 



sample with said spectroscopic device; 
(c) calculating a predicted set of analyte levels for each of 
the calibration algorithms in response to the set of 
non-invasive measurements, each of the predicted 
10 sets of analyte levels being characterized by a slope, a 

variability range, and a standard error of prediction; 
and 

ill (d) selecting an appropriate calibration algorithm based 

I f|- on the slope, the variability range, and the standard 

tf^ 15 error of prediction for each of the predicted sets of 

i^-; analyte levels. 

m 

5 

2. 'A method according to claim 1 wherein step (d) comprises 

tJl (i) selecting the predicted sets of analyte levels in wliich 

the standard error of prediction is less than an upper 
Sk3? 20 error limit, the variability range is greater than a 

lower range limit, and the slope is between a first 
lower slope limit which is less than 1 and a first 
upper slope limit which is greater than 1; 
(ii) for each of the predicted sets selected in step (i), 
25 calculating a suitability score in response to the slope 

and the standard error of prediction for that predicted 
set, and selecting the calibration algorithm 
corresponding to the predicted set having the optimal 
suitability score as the appropriate calibration 
30 algorithm; 
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(Ui) if no predicted ects arc selected in step (i), selecting 
the predicted sets of analytc levels in which the 
standard error of prediction is less than the upper 
error limit, the variability range is greater than the 
lower range limit, and the slope is between a second 
lower slope limit and a second upper slope limit 
which is less than or equal to said first lower slope 
limit; 

(iv) for each of the predicted sets selected in step (iii), 
calculating a suitability score in response to the slope 
and the standard error of prediction for that predicted 
set/ and selecting the calibration algorithm 
corresponding to the predicted set having the optimal 
suitability score as the appropriate calibration 
algorithm; 

(v) if no predicted sets are selected in step (i) or step (iii), 
selecting the predicted sets of analyle levels in which 
the variability range is lower than the lower range 
limit and in which the, standard error of prediction is 
less than the upper error limit/ 

(vi) from each of the predicted sets selected in step (v), 
selecting the calibration algorithm corresponding to 
the predicted set having the lowest standard error of 
prediction as the appropriate calibration algorithm; 
and 

(vii) if no predicted sete arc selected in step (i), step (in), or 
Step (v), determining that no calibration algorithm is 
appropriate. 

A method according to claim 1 wherein step (a) comprises 
(i) compiling non-invasive and corresponding reference 
data sets of measurements of analyte levels for each 
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of a number of samples; 

(ii) disregarding data sets thai are not suitable for 
calibration; 

(iii) combining data sets that are suitable for calibration 
into a plurality of groups depending on whether 
correlations of the combined data sets meet 
predetermined criteria; and 

(iv) generating o calibration algorithm for each of the 
groups? of data sets. 

4. A method according to claim 3 wherein the criteria in step (3) 
include minimizing or lowering correlations of the combined data sets In 
a particular group with parameters other than said analyte and 
maximizing or increasing the correlation between data sets in a particular 
group. 

5. A method according to claim 3 wherein steps (iii) and (iv) are 
performed using partial least-squares regression analysis. 

6. A method according to claim 1 wherein steps (a), (c), and (d) 
are performed on a computer associated with the spectroscopic device, 

7. A method according to chiim 1 wherein the analyte is glucosC/ 
the sample is an individual patient, and the spectroscopic device is a near- 
infrared spectrophotometer. 
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